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We investigate the nature of the ground state of the one-dimensional t — J model coupled to 
adiabatic phonons by use of the Lanczos technique at quarter filling. Due to the interplay be- 
tween electron-electron and electron-phonon interactions, the model undergoes instabilities toward 
the formation of lattice and charge modulations. Moderate on-site and intra-site electron-phonon 
couplings lead to a competition of different spin-Peierls and dimerized states. In the former case two 
electrons belong to the unit cell and we expect a paramagnetic band insulator state, while lattice 
dimerization leads to a Mott insulating state with quasi long range antiferromagnetic order. The 
zero temperature phase diagram is obtained as a function of intra-site and inter-site electron-phonon 
couplings, analytically in the J — > limit and numerically at finite J/t. 

PACS numbers: 75.10.Jm 05.70.Jk 75.40.Cx 



I. INTRODUCTION 

The coupling between electrons and lattice degrees 
of freedom in strongly correlated systems has been the 
subject of several experimental and theoretical inves- 
tigations showing that low dimensional materials are 
prone to structural distortions driven by electron-phonon 
interaction^^. The most celebrated case is the Peierls in- 
stability, extensively studied in weakly correlated mod- 
els: the ground state of a one dimensional electron- 
phonon system is characterized by a lattice distortion 
which opens a gap at the Fermi energy and gives rise 
to a charge modulation. In a half filled band, the den- 
sity wave is commensurate to the lattice, the unit cell 
doubles and the chain shows a spontaneous dimerized 
ground state. When on site Coulomb repulsion cannot 
be neglected, the system displays competition between a 
Mott insulator state with strong antiferromagnetic cor- 
relations, favored by electron-electron interaction, and a 
paramagnetic band insulator state stabilized by lattice 
effects. 

Experimentally, there are indications for the impor- 
tance of both electron-electron and electron-phonon in- 
teractions in quasi-one dimensional systems. In partic- 
ular there is a wide class of low dimensional materials, 
like the Bechgaard-Fabre salts series, where the presence 
of these interactions leads to several modulated phases 
such as: charge-density- wave (CDW), spin-Peierls (SP), 
antiferromagnetic (AF), spin-density- wave (SDW), and 
even a high temperature superconducting stated. 

The physical mechanism which drives these different 
orderings is not fully understood yet. A satisfactory the- 
ory should be able to explain the origin and the stability 
of different patterns, and the temperature scale at which 
they occur. In any case, the role of electron-phonon and 
electron-electron interactions will be instrumental for the 
explanation of the presence of bond and charge order. 

By now, the physics of one dimensional (ID) electron 
systems in a rigid lattice is well understood^ in terms of 
the Luttinger Liquid (LL) paradigm: a LL is a ID para- 
magnetic metal characterized by gapless excitations both 



in the spin and in the charge channel. In spin isotropic 
models, the low energy properties are determined by a 
single parameter denoted by K p , in particular K p < 1 
for repulsive interactions and K p > 1 in the attractive 
regime. The former case is characterized by strong anti- 
ferromagnetic correlations, while in the latter dominant 
superconducting fluctuations are present. Luttinger Liq- 
uid theory is an invaluable tool for the study of the ef- 
fects of perturbations on an interacting ID electron gas. 
Within this theoretical framework, it has been shown 
that electron-lattice coupling, impurities and interchain 
interaction can destabilize the LL behavior opening a gap 
in the spin and/or charge excitation spectrum. A charge 
gap leads to an insulating state whose properties depend 
whether the spin sector remains gapless or not. A gapless 
spin spectrum leads to a ID Mott insulator with (quasi) 
antiferromagnetic ordering, while a gap in the spin chan- 
nel corresponds to a band insulator. In short range mod- 
els, a spin gap with no charge gap preludes to a (quasi) 
superconducting state. 



Lattice effects in quasi one dimensional organic con- 
ductors have been mostly investigated on the basis of 
the Hubbard model with additional phonon couplings 
in adiabatic approximation (Peierls-Hubbard model)^. 
Recent works extended the study to include quantum ef- 
fects for the phonons^. However, the role of phonon cou- 
pling in ID models which may display attractive inter- 
actions, like the t— J model, received comparatively less 
attention^. The simultaneous presence of both bond 
and on-site (Holstcin) phonons is likely to be an essential 
ingredient for a faithful representation of quasi ID or- 
ganic materials. In this paper we explore systematically 
the effect of adding electron-phonon interactions to the 
t — J model by including inter-site (bond) and on-site 
(Holstcin) phonon couplings in the adiabatic approxima- 
tion. We aim at the understanding of the qualitative 
features of the lattice periodicities which are stabilized 
and the sequence of modulations the model undergoes 
when the strength of phonon couplings is varied. 
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II. THE MODEL 

The t — J hamiltonian provides one of the simplest 
models able to capture the essential low energy physics 
of doped antifcrromagncts. The electron dynamics is gov- 
erned by a competition between two terms: the kinetic 
contribution which describes hopping between adjacent 
sites and the nearest neighbor Hciscnbcrg interaction, 
which favors antifcrromagnctic alignment of spins. An 
infinite on site repulsion is also included in the model by 
means of a constraint. The phase diagram of the pure 
t — J model has been numerically determined^ while an- 
alytic solutions arc known only in the J — > limitii and 
at the supersymmetric point J/t = 3i2*i^. At half fill- 
ing the model reduces to the Heisenberg chain and when 
it is adiabatically coupled to the lattice via an arbitrar- 
ily weak perturbation, the system undergoes spontaneous 
dimerization (spin-Peierls instability) leading to a lattice 
distortion with wave vector q = 2k f'- adjacent spins are 
paired in a singlet state and a spin gap in the excitation 
spectrum shows up. For electron density n < 1 the pure 
t — J model enters the Luttinger Liquid phase and, at 
larger J/t, phase separation sets iniS. 

The one-dimensional t — J model coupled with adia- 
batic phonons is defined by the hamiltonian 

H = -^(l-* i )(4« r C^l,a + /».C.) + 
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where Si are spin-1/2 operators at the site z, cj — 
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are electron Gutzwiller-projected cre- 



ation operators and = c\ a Ci t(T is the local electron 
density. We have also set the bare hopping integral t to 
unity thereby fixing the energy scale. Hamiltonian 
includes both nearest neighbor superexchange coupling 
J and Coulomb repulsion V. The hamiltonian paramet- 
rically depends on the classical variables Si and i>j which 
identify the bond distortion and the amplitude of the in- 
ternal site vibration, respectively. The bond distortions 
arc defined as Si = (v-i — itj_i), where m is the displace- 
ment of the i site from the equilibrium position, while 
the on-site displacement u, corresponds to an effective 
parameter which combines several processes. Clearly, the 
physical range where the hamiltonian Q may be used is 
restricted to |<$$| <C 1. The first three terms correspond to 
the usual t — J — V model coupled with adiabatic bond 
phonons, the next one is the Holstein type interaction 
between the particles and the lattice, while the last two 
terms describe the elastic deformation energy, Kb and 
Kh being the spring constants. 



The t — J model can be regarded as the strong cou- 
pling limit of the extended Hubbard model when U — > oo. 
This leads to a definite relationship between the Hubbard 
and t — J couplings: J = 4/ (U — V) and g = 2&. However, 
it is knownia that the t — J model can be also consid- 
ered as an effective one band hamiltonian unrelated to the 
Hubbard model. In this case, the superexchange coupling 
J takes into account electronic processes involving other 
degrees of freedom (like the Oxygen orbitals in Copper 
Oxide materials and in manganitcsJ^) . The study of the 
effects of phoiion coupling in the t — J model allows to 
analyze the role of magnetic interactions in stabilizing 
lattice distortions in itinerant electron models. By vary- 
ing the J parameter we can explicitly tune the strength of 
antifcrromagnctic coupling in the system: the compari- 
son between the (analytical) J — > limit and the (numer- 
ical) finite J results will contribute to the understanding 
of the physical mechanisms leading to the observed pe- 
riodicities. In the following we will therefore consider J, 
V and g as free parameters in the t — J hamiltonian • 

The calculations have been performed at quarter fill- 
ing n = 1/2, which is an appropriate choice to mimic the 
behavior of the organic compounds. By use of Lanczos 
diagonalizations, we found the lattice distortions provid- 
ing the lowest energy, without any assumption on the 
periodicity of the broken symmetry ground state. The 
extremum condition defining the optimal modulation can 
be formally expressed as 
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The sum rule J2i $i — reflecting the ring geometry of 
the model has been implemented by use of a Lagrange 
multiplier, while no constraint on {vi} has been imposed. 
Here (...) is the ground-state expectation value obtained 
by exact diagonalization using the Lanczos algorithm. In 
order to find the bond length configuration which mini- 
mizes the total energy E({Si} , {fj}), at fixed spring con- 
stants (Kb, Kh), we iterate the set of equations J2J) until 
convergence is reached. 

Following a previous investigation^ we chose open shell 
boundary conditions in order to mimic the gapless behav- 
ior of the (pure) t — J model in the thermodynamic limit 
when studying a small size chain (N < 16). 



III. ANALYTICAL RESULTS 

In the special limit J — > (and V = 0) the problem can 
be solved analytically even in the thermodynamic limit 
because the charge and spin degrees of freedom exactly 
factorize at all length scales^. By setting J = in Eq. 
Q the hamiltonian reduces to the pure kinetic contribu- 
tion. In ID the single occupancy constraint of the t — J 
model can be explicitly taken into account by promot- 
ing the spinful electron to spinless fermions. The origi- 
nal electron problem at quarter filling then maps into the 
adiabatic Su-Schricffcr-Hccgcr (SSH) modcU^ for spinless 
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fermions at half filling with an additional local potential 
Vi and elastic deformation energies. 

Being at half filling, the hamiltonian 0) is unstable 
toward dimerization and a gap in the energy spectrum 
opens for arbitrarily weak electron phonon coupling. In 
the SSH adiabatic case («, = 0) a dimerized lattice distor- 
tion pattern with uniform charge is stabilized and gives 
rise to a bond ordered wave (BOW) . In the pure Holstein 
case (Si = 0) the lattice is undistorted and a CDW forms. 
A competition between BOW and CDW dimerizations is 
expected in the general case. By assuming a dimerization 
pattern of the form Si = (— l) l S and Vi = (— l) l v — vo> 
the hamiltonian for J — ► becomes: 

H = - J2 (1 - *("!)*) ('•;..' -• i - + h.c) + 

i 

+ v J2(-iyni+^K B 5 2 + 

i 

+ yA'jyi; 2 - v N e + ^K H v% (3) 

where N is the number of sites and N e — N/2 is the 
number of particles. 

In this case the hamiltonian can be easily diagonalized 
by canonical transformations and the total energy, in the 
thermodynamic limit, has the form: 

E = —N J*'" |j- [4 (cos 2 q + S 2 sin 2 q) + v 2 } h + 

The minima of the energy A as a function of the dimer- 
ization parameters 5 and v always correspond to pure 
BOW (v = 0) or pure CDW (S = 0) according to the 
values of the elastic constants. The two modulations are 
respectively given by the equations: 

f 77 / 2 fjn _1 

K B = / — 2 sin 2 q (cos 2 q + S 2 sin 2 q) 2 
K H = ^-(Acos 2 q + v 2 )^ (5) 

J-ir/2 Zn 

Finally, the BOW/CDW phase boundary is obtained by 
equating the energies of the two distortions. The result 
is shown in Fig^ the asymptotic form of the transition 
line for large elastic constants (i.e. for weak distortions) is 
given analytically by Kb = 4Kb — 2/ir. We observe that 
at small Kb a charge dimerization in the undistorted 
lattice (Ak p CDW state: see FigEJl) prevails, while in 
the rest of the phase diagram a dimerized lattice with 
uniform charge (4k p BOW state: see Fig|2t) is stabilized. 
In both cases a gap in the charge excitation spectrum is 
present and the system is therefore an insulator. 

In order to characterize the electronic ground state of 
the model we observe that even in the presence of adia- 
batic phonons the spin-charge factorizationii holds and 




l/K, 



FIG. 1: Zero temperature phase diagram of the t — J model 
coupled to adiabatic phonons in the J — » limit. 

therefore, following the analysis carried out for the strong 
coupling limit of the one dimensional Hubbard modelii, 
the spin correlations can be expressed as a convolution^: 

■r + l 

(S r -S ) = Y^P r SFij)SH(j-l) (6) 

where Pg F (j) is the probability of finding j particles in 
(0, r) with one particle in and another in r, evaluated 
in the ground state of the dimerized spinlcss Fermi gas. 
If we set iV r = X)[=o n i: tnc probability Pg F (j) can be 
formally expressed as Pg F (j) = (n n r S (N r — j)). 

The asymptotic decay of Eq. @), in the presence of 
lattice dimerization and uniform charge (Fig[3Jt) has been 
previously evaluated in Ref A: 

cos(2k F r -?) 
Sr • S ) oc — { — *L (7) 

r 

and shows antifcrromagnetic quasi long range order. 
Here the Fermi wave vector is given by kp = nn/2 = tt/4. 
An analogous scenario is expected in the CDW phase, 
where the lattice is uniform but the charge is dimerized 
(see Fig|2JI). In this case the spin-spin correlation func- 
tion is not translationally invariant and decays as: 

(S r • So) oc C -^^l (8) 

r 

within each sublatticc. This asymptotic result, valid for 
any charge dimerization, can be obtained in a simple 
way in the limiting case where the electron density is 
zero in one sublattice (i.e. when v — ► oo). In this case 
Psf (j) = <5j",r/2+i which leads to © via Eq. |JBJ. There- 
fore in both the BOW and CDW case, the t — J model 
is characterized by a charge gap together with a 1/r de- 
cay of spin correlations and diverging antifcrromagnetic 
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FIG. 2: Pictorial representations of the competing density 
waves in the ground state of the t—J model, a) 4&f BOW two 
different bonds and uniform charge, b) 2k f BCDW tetramer- 
ized phase (I-S-I-W) three different bonds and two different 
charges, c) 4k F CDW-SP (W-W-S-S) two different bonds and 
three different charges, d) 4&f CDW uniform bond and two 
different charges. 

susceptibility, as in the ID Hciscnbcrg model. Such a 
scenario strongly suggests that, even for J — > 0, lattice 
coupling favors antiferromagnetism, which will be even- 
tually stabilized when the inter-chain coupling, present 
in the three dimensional systems, is taken into account. 

Interestingly, this solution directly applies to the 
strong coupling limit of the Hubbard model via the corre- 
spondence J = 4/U. Therefore we expect that the phase 
diagram of Fig. ^faithfully represents the behavior of the 
U — > oo Hubbard model when only the on site Coulomb 
repulsion is considered (i.e. V = 0). The presence of fi- 
nite range interactions is likely to play an important role 
in stabilizing other periodicities. 

IV. LANCZOS DIAGONALIZATIONS 

We now turn to the analysis of the full hamiltonian (|TJ 
including the super-exchange term J and both Holstein 
and Peierls lattice couplings. Unfortunately for finite val- 
ues of J it is not possible to derive exact results in the 
thermodynamic limit and then we investigate this model 
at quarter filling and vanishing magnetization. Results 
are given for a fixed value of magnetic and magnetoelastic 
couplings (J = 1.4, V = and g = 1) in order to dis- 
cuss the interplay between bond and charge distortions in 
the presence of antifcrromagnctic interactions. We note 
that, due to strong finite size effects, a large value of the 
magnetic coupling J is necessary to stabilize the possible 
symmetry broken phases. Otherwise, lattice and charge 
modulations are present only for exceedingly small elas- 
tic constants. However, we believe that the qualitative 
features of the phase diagram at realistic values of J are 
well reproduced by our small cluster calculations, while 
we expect that the position of the phase boundaries and 
the amplitude of the lattice distortions will be severely 
renormalized in the thermodynamic limit. 
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FIG. 3: Phase diagram of the t-J model coupled to adi- 
abatic phonons Q obtained by numerical diagonalization 
of a N < 16 site chain. The boundaries of the dimerized 
bond order 4&f BOW, the tetramerized 2&f BCDW, the 4&f 
CDW and uniform bond order, the 4&f CDW-SP tetramer- 
ized bond-charge order are schematically shown. The dashed 
line represents a crossover between two different metastable 
states, as discussed in the text. 



In the ground state we observe only two competing pe- 
riodicities: tetramerization (g = 2kp) and dimerization 
(q = 4kp). In the former case two electrons per unit 
cell are present and then the system will probably be a 
paramagnetic band insulator with both charge and spin 
gap. This regime corresponds to the spin-Peierls phase 
experimentally found in the Fabre and Bechgaard salts 
series^. When a q = Akp modulation is stabilized, two 
sites belong to the unit cell thereby accommodating a 
single electron: the Umklapp scattering opens a charge 
gap and the system behaves like a Mott insulator. 

In FigO we present a schematic ground state phase 
diagram in the 1/Kh-1/Kb plane, showing the regions 
where the various broken symmetries dominate. Note 
that, as soon as both phonon couplings are significant, 
mixed CDW-BOW ground-states are stabilized. 

All numerical results have been obtained by Lanc- 
zos diagonalization with no assumption on the super-cell 
structure or on the periodicity of the ground state. This 
limits the investigated cluster sizes to the range N < 16. 
In all cases the results do not depend on the starting 
bond configuration nor on the convergence requirement 
of the adopted iterative method^. At large Kh, the nu- 
merical data show a narrow region where the ground 
state is a 4/cp BOW, as we found in the J — > limit. 
However, we suspect that in the thermodynamic limit 
this phase disappears at finite J: our evidence for this 
comes from few Lanczos diagonalizations carried out in a 
N = 20 site chaini^ showing that the stability region of 
the Akp BOW phase is considerably reduced, while the 
other phase boundaries hardly move. Therefore we can 
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conclude, as already suggested for the Peierls-extended 
Hubbard models, that the Akp BOW ground state phase 
could be an artifact of the limited size chain. 

Few calculations have been also performed for differ- 
ent parameter values. In all cases we observed the same 
sequence of states shown in the phase diagram of Fig. 
[3] while only the phase boundaries are sensibly affected 
by the precise values of the couplings: When J is re- 
duced, a uniform phase appears at large values of Kh 
and the stability domains of the BOW phase gets larger. 
Remarkably, the qualitative appearance of the phase di- 
agram remains unaltered by tuning the magnctoelastic 
coupling g and most of the changes can be absorbed into 
a renormalization of the couplings, the effective J being 
reduced when g is increased. Similarly, by including the 
nearest neighbor Coulomb repulsion V > the rigidity of 
the system is enhanced and the effective elastic constants 
get larger. However, no qualitatively different behaviors 
appear. 

For the investigated value of the super-exchange term 
J, all the various charge and bond modulations can be 
stabilized by tuning the Holstcin coupling at fixed Kb- 
Moving from weak to strong on-site coupling (i.e. by low- 
ering Kh) three different types of structures are found: 
a Akp CDW, a Akp CDW-SP and a 2k F BCDW. In all 
these phases, the bond and charge density pattern can 
be simply parameterized as: 



BCDW | 4K r CDW-SP 



4K, CDW 



{rii) 



1 



Aw.*- cos (iri) + A 2 k^ cos —i 



Si = B ikp cos (tt£) + B 2kp cos ( —i + <fo 



(9) 



At small on-site coupling (Kh — > oo) Akp CDW instabil- 
ity prevails, while bonds remain uniform (see Fig|2Ji). In 
this phase, the only non vanishing amplitude in Eq. @ 
is Aik F ■ Note that a previous investigation of the t — J 
model without Holstein phonons^ showed no evidence in 
favor of the Akp CDW-SP and Akp CDW phases. By low- 
ering Kjj, a tetramerized phase denoted Akp CDW-SP is 
obtained in the region of intermediate phonon couplings 
(see FigEt): the charge ordering is characterized by a 
superposition of a 2kp and a Akp density wave, while 
the lattice modulation contains only a strong 2k p com- 
ponent. This regime corresponds to the W-W-S-S bond 
sequence, where S and W respectively represent a weak 
(W) and strong (S) bond. The charges order according 
to the sequence 1-0-1-0, where 1 (0) identifies a local den- 
sities higher (lower) than the average (n = 1/2). Note 
that, as shown in FigEfc, the two low density ("0") sites 
inside the unit cell have different densities. The bond 
modulation is reproduced by setting B 2 k F ^ 0, B^k F — 
and (j>b — tt/4 in Eq. @, while the charge pattern re- 
quires A 2 k F ^ 0, A±k F ^ and C = 0. At small Kh the 
ground state is characterized by a mixed 2kp and Akp 
lattice order accompanied by a strong 2kp charge modu- 
lation (BCDW phase). The bond sequence is of the form 
I-S-I-W (where I stands for intermediate, i.e. a mod- 
erately strong bond) and the charge order is the type of 
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FIG. 4: Local densities in the unit cell as a function of Kb 
at the given value Kh = 1 



0-1-1-0 (see Fig|2b). The site centered BOW corresponds 
to both Bik F ^ and B 2 k F ^ with fa = 0, and the 
2kp bond centered CDW is parameterized by A 2 k F ^ 0, 
A^k F — and <j> c = tt/4. Notice that while the strong (S) 
and weak (W) bonds are affected by the strength of the 
on-site and bond couplings, the intermediate (I) bonds 
remain roughly constant throughout the BCDW region. 
By comparing these patterns to those obtained in the 
absence of Holstein phonons^ we see that in the BCDW 
phase the S-W-S-W bond sequence is never stabilized 
(except in a narrow region close to the phase bound- 
ary). However, by reducing the antifcrromagnctic cou- 
pling or by including a nearest neighbor Coulomb term 
V, a small BCDW region characterized by the S-W-S- 
W sequence does appear, as already noticed in a pre- 
vious investigation^. This suggests that in the small J 
limit, where the t — J model reduces to the strong cou- 
pling regime of the Hubbard model, the S-W-S-W bond 
pattern shows up^ 



The sequence of modulations displayed in the phase 
diagram of Fig|21does not change qualitatively with Kb 
and is maintained even at small bond phonon coupling, 
i.e. in the limit Kb — > oo, where the amplitude of the 
bond distortions becomes vanishingly small (Si — > 0). 
Contrary to the Hubbard model, the quantitative val- 
ues of the site charge density wave (n,) in the Akp 
CDW-SP and Akp CDW state are quite small, whereas 
in the BCDW state they are considerably larger. We 
could not stabilize, in our calculations, any Akp CDW- 
SP or Akp CDW states with large values of (rii), nor any 
BCDW state with small (m). In FigH we show a typical 
trend of the site charge density wave {n%). The vari- 
ous points at each Kb provide the independent density 
values inside the unit cell of the superlattice generated 
by the phonon distortions. Both the Akp CDW-SP and 
the Akp CDW regime are characterized by a substan- 
tially i^s-independent value of the highest density. In- 
stead, the lower density of the Akp CDW state splits into 
two branches in the Akp CDW-SP phase: the splitting 
increases continuously even through the BCDW phase 
boundary. Similarly, by varying Kh at fixed Kb the size 
of the bond modulation and its periodicities are shown 
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FIG. 5: Bond strengths in the unit cell as a function of Kh 
for K B = 3 

in Fig. [3] in the different regions of the phase diagram. 

Whenever lattice tctramcrization is stabilized, two 
electrons belong to the unit cell and then the chain be- 
haves as a paramagnetic band insulator, like in the spin- 
Peierls phase. We believe that, because of the presence of 
both charge and spin gap, this phase is stable toward 3D 
interchain coupling. While for lattice dimcrization one 
electron belongs to the unit cell: electron repulsive inter- 
actions drive the system toward a Mott insulator phase 
and, if 3D coupling is allowed, antiferromagnetic order is 
favored. 

Our numerical study has been limited to investigate 
the ground state properties of the hamiltonian Q . How- 
ever a finite temperature analysis is relevant in order to 
apply the results of our model to real materials. Some 
preliminary information about the temperature depen- 
dence can be inferred by evaluating, via Lanczos diago- 
nalizations, the energies of the lowest metastable state of 
the hamiltonian Q obtained by imposing a given period- 
icity of the bond/charge modulation. More precisely, let 
us consider a point of the phase diagram corresponding 
to a 2kp modulation of the ground state, either of the 
BCDW or of the 4fc F CDW-SP type: We performed a 
self consistent Lanczos diagonalization by imposing the 
further constraint of Akp periodicity, i.e. of a two site 
unitary cell. In this way, the calculation cannot converge 
to the actual ground state (which is characterized by a 
four site unitary cell) but instead it selects a metastable 
state identified as the lowest energy state with two site 
periodicity. The nature of this metastable state changes 
throughout the phase diagram of Fig. [3J it is of the Ak F 
CDW type (Fig. Hi) in the whole Ak F CDW-SP region 
and to the right of the dashed line in the BCDW phase, 
while it is of the Ak F BOW type (Fig. in the re- 
maining part of the BCDW region. The picture of the 
underlying energy landscape and the evaluation of the en- 
ergy barrier separating the symmetry broken states can 
be also obtained by Lanczos diagonalizations. Starting 
from a W-W-S-S configuration (Fig. H) we gradually 
changed the difference 5s — 5w between the two strong 
and the two weak lattice distortions in the unitary cell, 
while letting the self consistent Lanczos algorithm opti- 
mize the vibrational parameters Vi. A typical plot of the 
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FIG. 6: Energy landscapes obtained by varying the bond 
length difference 5s — Sw (see text) for Kb = 2, Kh = 2 
(upper panel) and Kb = 5, Kh = 1 (lower panel). 



energy as a function of the bond asymmetry is shown 
in Fig. H3 (lower panel). The broken symmetry Ak F 
CDW-SP ground states correspond to the energy min- 
ima at non vanishing bond asymmetry, while the 4k p 
CDW metastable state with uniform bonds is the central 
maximum. When the thermal energy ksT is higher than 
the energy barrier AE we expect that thermal fluctua- 
tions will restore the full symmetry among the bonds in 
the unitary cell stabilizing a Ak F CDW state. Analogous 
calculation can be performed in the other regimes. In 
Fig. El (upper panel) the energy landscape corresponding 
to the BCDW phase I-S-I-W is shown. Here, 5s — 5\y 
identifies the difference between the strongest (S) and 
the weakest (W) bond, while the two other bonds (of 
equal strength) are kept constant. The overall picture is 
quite similar to the previous case, although the energy 
barrier is considerably higher. The metastable state at 
$s = $w now corresponds to a dimerized (BOW) configu- 
ration with uniform charge: such a phase is then expected 
to represent the intermediate temperature regime in this 
region of the phase diagram. At higher temperatures a 
metallic uniform state is stabilized. The two scenarios 
discussed here are schematically shown in Fig. [7| In- 
terestingly, the former case (panel a) corresponds to the 
ground state sequence experimentally found in the Fabrc 
salt series by increasing the temperature^. 



V. CONCLUSIONS 

We have investigated the t — J hamiltonian includ- 
ing both adiabatic Holstein and Pcicrls lattice couplings, 
with amplitudes defined by the two independent param- 
eters 1/Kjj and 1/Kg, respectively. Exact analytical 
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FIG. 7: Schematic expected sequences of phases obtained 
from a SP ground state (either BCDW or 4k F CDW-SP) by 
increasing the temperature 



treatment of the J — ► limit together with the im- 
plementation of a numerical method based on Lanczos 
technique allowed to determine the ground state phase 
diagram as a function of the elastic constants at fixed 
value of the superexchange interaction J = 1.4. For a 
weak on-site (Holstein) coupling (i.e. Kh — * oo), two 
different dimerized regimes are stabilized: a 4k f BOW 
and a 4fe CDW which display the typical features of 
a Mott insulator and suggest antiferromagnetic order in 
three dimensional samples. In the former case the charge 
distribution remains uniform while the lattice is dimer- 
ized; in the latter case the bond is uniform with dimerized 
charge. When Kh gets smaller, tetramerized phases ap- 
pear displaying the typical behavior of the spin-Peierls 
regime: the 4k,F CDW-SP phase, for intermediate cou- 
pling, and the Ikp BCDW for strong coupling. The SP 
phase with 1-0-1-0 charge order has two kinds of bond 
and three distinct charges, the bond distortion pattern 
is W-W-S-S and the distortion makes the charges on the 
sites labeled "0" unequal (4kp CDW-SP). In contrast, 
the SP phase with 0-1-1-0 charge order is a BCDW phase, 
with two different charges and three different bond order, 
the bond distortion pattern is I-S-I-W. In the tetramer- 
ized phases two electrons belong to the unit cell and then 
we expect a paramagnetic band insulator, and we be- 
lieve that because of the presence of both charge and 
spin gap, these two SP phases are stable toward 3D cou- 
pling and faithfully represent the low temperature phases 



of the model. The SP regimes give rise to two differ- 
ent scenarios at higher temperature: by increasing the 
temperature, the system will explore metastable config- 
urations corresponding to different bond/charge pattern 
before entering a paramagnetic undistorted phase. By 
examining the metastable states of our hamiltonian we 
found two possibilities, one of them, the sequence showed 
in Fig[7^, corresponds to the experimental phase diagram 
of Fabre salts in the temperature-pressure planed. Our 
study predicts that another sequence of phases, shown 
in FigUJa, should be possible in quasi one dimensional 
materials characterized by larger Kh- 

The effects of charge and bond ordering have 
been extensively investigated in the extended Hubbard 
modeii2i3i adiabatically coupled with the lattice. A 
ground state phase diagram obtained as a function of 
the electron-phonon interaction shows the same period- 
icities that we found in our model, with the exception 
of the BCDW regime whose bond pattern follows the 
S-W-S-W' sequence. By lowering J or by adding the 
nearest neighbor Coulomb repulsion V, the S-W-S-W' is 
recovered in a small window close to the BCDW phase 
boundary. 

This analytical and numerical study represents the first 
attempt to investigate the interplay of bond and Holstein 
phonons in the one dimensional t — J model. The re- 
sults show an overall similarity with respect to previous 
studies for the Hubbard model. However, even when the 
bond pattern is similar, the corresponding charge density 
wave displays different features between the two models 
due to the presence of a finite range attraction, of mag- 
netic origin, in the t — J systems. A nearest neighbor 
attraction smears the strong charge oscillations present, 
for instance, in the 4k f CDW phase of the Hubbard 
models. A direct measure of charge and bond patterns 
in organic conductors will help in understanding the role 
of direct magnetic coupling in real materials. Finally, 
another remarkable result which emerges form this work 
is the counter-intuitive suppression of antiferromagnetic 
(quasi) long range order in strongly correlated models 
when J is increased, due to the coupling to the lattice 
which drives a 2kp instability, leading to a Spin-Peiers 
state. 

We acknowledge fruitful discussions with V. Lante. 
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PRIN grant. 
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